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SUMMARY 

A previous analysis of turbulent flow and heat transfer in smooth 
tubes -irtth variable fluid properties was modified in order to be" appli- 
cable to air -which heis a Prandtl mmiber slightly less than unity (0.73). 
In order to verify the analysis^ tests were conducted to determine 
local heat -transfer coefficients and friction factors for fully developed 
tixrbulent flow of air in a smooth electrically heated tube having an 
inside diameter of 0.87 inch and a length of 87 inches. The tests were 
conducted at high ratios of wall to fluid bulk temperature. Velocity 
and tenperature profiles were measured for some conditions. 

The analytical ly predicted results were found to agree closely -with 
the experimental. Bo-th the analytical and experimental resxilts indi- 
cated that the effects of ratio of wall to bulk temperature on the 
Nusselt number correlation can be eliminated by evaluating the fluid 
properties, including density, in the Reynolds and Nusselt numbers at 
a temperature close to the average of the wall and bulk tenperatures . 

In order for this correlation to apply it is necessary to assiime con- 
stant Prandtl number and constant specific heat in calculating both the 
experimental and analytical results. 

The effects of variation of shear stress and heat transfer across ' 
the tube on the velocity and temperature distributions were analytically 
investigated and found to be small for turbulent flow -with variable 
fluid properties. The effects of molecular shear stress and heat trans- 
fer in the region at a distance from the wall were also investigated 
and found to be small. 


' INTRODUCTION 

In an experimental investigation of average heat transfer and 
friction coefficients for air flo-wing in smooth tubes (reference l), it 
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was found that the ratio of wall to fliiid hulh temperature has an appre- 
ciable effect on the heat -transfer and friction correlations . A theo- 
retical analysis given in reference 2 indicates a similar effect for local 
coefficients. For verification of the analysis given in reference 2 , it 
is desirable to obtain local rather than average heat-transfer coeffi- 
cients and friction factors and meastarements of radial velocity and tem- 
perature distributions for large heat-transfer rates . Local rather than 
average coefficients are desirable because average coefficients include 
possible effects of flow development near the tube entrance and of heat 
loss at the enfla of the tube. 

In the investigation reported herein^ which was conducted at the 
KACA Lewis laboratory^ local heat-transfer coefficients and friction 
factors at hl^ ratios of w all to buli temperature were measured at a 
point in a tube I'diere the flow was fully developed. Velocity and tem- 
perature distributions were also measured for some conditions . The 
analysis for flow and heat transfer with variable fluid properties given 
in reference 2 is modified for application to air, which has a Prandtl 
number sli^tly less than 1, and the analytically predicted results are 
compared with the data. 


ANALYSIS FOE PEANDTL NUMBERS DIFFERING SLIGHTLY FROM ONE 


QTe analysis given in reference 2 is for a fltild with a Prandtl 
number of 1. Inasmuch as air has a Prandtl number of about 0.73, the 
analysis is modified in this section in order to make a comparison 
between analytical and experimental results. 

For obtaining the velocities and temperatures in the tube as func- 
tions of distance from the wall, the differential equations for shear 
stress and heat transfer are often written in the following form: 

du du 

t: = IJ- — t pt w— 

^ dy dy 



, dt c. 

<1 = - ^ - P V'n ^ 

(The symbols used in this report are defined in appendix A.) 
The second equation can be rewritten as 
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In the preceding eq.nations, 6 and are the coefficients of eddy dif 

fusivity for momentum and heat transfer, respectively, the values for 
■which are dependent on the amount and kind of turbulent mixing at a 
point. 

Assumptions . - 0!he following assumptions are made in the use of 
equations (l) and (2) for obtaining velocity and temperature distri- 
butions with heat transfer; 

(1) The eddy diffusivities for momentimi and heat transfer e and 
are equal. Previous analyses for flow in tubes based on this assump 

tion yielded heat-transfer coefficients and friction factors that agree 
I'rLth experiment. 

( 2 ) The expressions for eddy diffusivity that are found in refer- 
ence 3 to apply to flow without heat transfer apply also to flow with 
heat transfer -S'rith variable fluid properties.- These expressions are 

e = n^uy (3) 

for flow close to the wall (y-4-< 26} and the Karm^ relation 



for flow at a distance from the wall (y+>26), where n and K are 
experimental constants ha-ying the values 0.109 and 0.36, respectively. 

( 3 ) The variations across the tube of the shear stress t and the 
heat transfer per uni.t area q have a negligible effect on the velocity 
and temperature distributions. Ihe effect of these variations will be 
investigated in appendix B and in the section "Predicted Effect of Var- 
iation of Shear Stress and Heat Transfer Across Tube." 

( 4 ) The molecular shear stress and heat-transfer terms in equa- 
tions ( 1 ) and ( 2 } can be neglected in the region at a distance from the 
wall. The effect of these terms will be investigated in appendix C and 
in the section "Predicted Effect of Mcleoilar Shear Stress and Heat 
Transfer in Region at a Distance from Wall." 

( 5 ) The static pressure can be considered constant across the tube. 

(6) The Prandtl number and specific heat can be considered constant 
with temperature variation. The variations id.th temperatxare of the spe- 
cific heat and the Prandtl number of gases are of a lower, order of mag- 
nitude than the variations of viscosity, thermal conductivity, and 
density. 
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Flov close to vaill . - For obtaining the velocity and temperature 
distrihutlons close to a smooth wall,' the e^^ression for e from equa- 
tion ( 3 ) is substituted into equations (l) and ( 2 ) to give 

Tq = (ti + n2puy) ^ 


and 


Cpg 



n^puy 


\dt 

7aSF 


idiere e and have been assumed eq.ua]. (assumption (l)) and r and 

q have been replaced by 'Tq and q^ (assumption (3)). With substitu- 
tion of the dimensionless quantities u"*", y"^, and p these equations 
become 


and 



(5) 


( 6 ) 


From viscosity data it is found that p/lig can be represented approxi- 
mately by (t/tQ)*^j and from the assumption of constant static pressure 
(assumption (5)) and the perfect gas law, p/pQ can be replaced by tQ/t. 
From the definitions of p and t**" 



( 9 ) 
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Equations (S) and (9) can "be -v/rltten in integral form as 


r\y 


+ 

u = 


dy 


0 -2^ 


1-pt 


and 


5 


( 10 ) 


t+ = 


ny 


lo 


Pr 


dy 

+Ad 


^ + + 
l-pf*” 


(U) 


Equations (lO) and (U) can le solved simultaneously "by iteration, 
that is, assumed values for u+, y"**, and t+ are substituted into the 
right sides of the equations and new values of u"*" and f*" are calcu- 
lated hy numerical integration. Olhese new values are then substituted 
into the ri^t sides of the equations and the process is repeated until 
the values of u"*" and f*" corresponding to each value of y"*” do not 
change appreciably. Equations (lO) and (n) give the relations between 
u"^, f*", and y+ for various values of the heat-transfer parameter p 
for flow close to a wall. 

Flow at a dist^ce from wall . - By making use of assumptions (l), 
(2), (3)j (4:), and ^6), eqiiations (l) and (2) can be written in dimen- 
sionless form for flow at a distance from the wall as 


1 = -£- (du~*~/dy+)^ du"*" 

PO (dV/dy+^)^ ^ 

and 

1 = -P. (du'*'/dy'*')^ dt'*' 

PQ (d^+/dy+2)2 dy+ 


( 12 ) 


(13) 


Dividing equation (12) by equation (l3) and integrating result in 




Au’’ 


dt+ = 


vjt^ 


du"* 


lui 
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or 

f*" - ti'*' = u+ - U1+ (14) 

t 

Tdiere nnfl are tlie values of and u"*" at ■wiiich is 

tile lowest value of for 'wtiicli equations (12) and (13) apply . Then 

= 1 - (u'^-U3^‘*'+tx^ P 

or^ from assumption (5) and the perfect gas law^ 


_P_ ^ 1 

PQ 1 - (u+-ux++t]^‘*’)p 


(16) 


On substituting equation (16) in equation (12) there, results 


1 


- (u^-Ux^+tx^) p = 


(du^dy"^)^ 

(d^uVdy"^^)^ 


One integration of equation (l7) gives 



(17) 


(18) 


where the negative sign was selected in taking the square root in order 
to make x positive. Integration of equation (l8) gives 


+ 

y = 


CxP 


2x 


^ a/i - p(u'*'-ux'^+tx'*') r 1 

^ ^ ^ - P(u-^-uiW) + ij 


+ c 


(19) 


As the wall is approached, the velocity gradient becomes very large com- 
pared with that at a distance from tbe wall so that dy+/du+ approaches 

zero' as y"*" approaches zero. If dy‘*'/du'*' is substituted from equa- 
tion (18) into equation (19) and dy+/du+ is set equal to zero idien y+ 
and u"*" equal zero, G equals zero. To determine Cx^ set u"*" 
when y+ = yx"^» Then 
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Equation (20) gives the relation "between u"*" and y"^ for various values 
of p for flow far from the wall. The parameter and thus the tem- 

perature distribution can he calculated from equation (14) . For a 

Prandtl number -of 1, t^^ = U3_'^j and equation (20) reduces to the equation 
for a Prandtl number of 1 given in reference 2. For p = 0 equa- 
tion (20) is indeterminate. For this case p is set equal to zero in 
equation (17) before integrating^ and the well-known logarithmic equa- 
tion is obtained; 


u+ = i loge y"^ + C 


Husselt number, Reynolds number, and friction factor . - It can 
easily be shown from the definitions of the quantities involved (see 
reference 2) that the Husselt number^ Reynolds nuoiber^ nnil friction fac- 
tor with the fluid properties evaluated at the wall temperature are 
given by 


Nuq 


V 


Rep 


2^bV 



( 21 ) 


where 


( 22 ) 

( 23 ) 
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(24) 


and 





- +2 

0 u 


(ro'^-y’^)dy'^ 


(25) 


The relations "between u"*", y'*', and f*" are calculated from the equa- 
tions given in the tiro preceding sections . 


The Nusselt numher^ Reynolds number, and friction factor with the 
fluid properties evaluated at some temperature in the fluid other than 
the wall temperature can be foimd by using the following equation for 





(26) 


Both the thermal-conductivity ratio and the viscosity ratio are equal 
to the temperature ratio raised to the exponent d inasmuch as the 
specific heat and Prandtl number are assumed constant (assumption (s)). 


Although the preceding analysis was carried out for the case i&ere 
the compressibility effects due to high velocities can be neglected, it 
is shoim in reference 2 that the same correlations apply, in general, 
to flow at high subsonic velocities. For the latter case (except for 
very small temperature differences) it is necessary only to replace the 
static bulk temperatvire by the total bulk temperature in the definition 
of the heat-transfer coefficient h. 
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APPAKATUS 

A schematic dieigram of the test section and associated equipment Is 
shovn In figure 1. Compressed air first passes through a filter and then 
Into an alumina- type air dryer. From the dryer the air flows through 
one of four orifices selected renotely hy orifice selection Talves and 
then through valves 'vdilch automatically control the flow rate. From 
the control valves the air flows Into a surge tahk, then to the Inlet 
mixing tank^ throvigh the test section, and Into the outlet mixing tank 
from which It Is discharged to the atmosphere through a pressure regulat- 
ing valve . 

Electric power Is supplied to the heater tuhe from a 208-volt 60-cycle 
svipply through, an autotransformer and a step- down transformer with 
maximum secondary voltages of either 15 or 25 volts . The low-voltage 
leads of the step-down transfoimer are connected to the flanges on the 
ends of the test section hy copper hus hars. The capacity of the elec- 
tric equipment Is 100 kilovolt-amperes. Additional heat to compensate 
for end losses Is supplied to each end of the heater tuhe -vdiere guard 
heater colls of nlchrome wire are wrapped around the cone-shaped heater- 
tuhe flanges. Individual control of the 110-volt 60-cycle supply to each 
end of the tuhe Is obtained hy variable transformer with capacities of 
1 kilovolt-ampere. 


Test Section 

The test section, shown In detail In figure 2, Is a commercially 
smooth Inconel tuhe having an Inside diameter of 0.87 Inch and an out- 
side diameter of 1,25 Inches -trLth a sharp right-angle entrance. Steel 
flanges with cone-shaped* centers welded to the tuhe at each end provide 
electrical contact with the transformer leads from the power supply. 

The tuhe, which Is thermally Insulated hy a 5-lnch thickness of high 
temperature quarts wool, acts as a heating element with an effective 
heat-transfer length of 87 Inches. Static pressures are measured 
through 0.03-lnch holes drilled In the ttibe at the positions shown In 
figure 2. Outside txibe-wall temperatures are measured hy chromel-alumel 
thermocouples (two thermocoiiples located 180*^ apart at each of 21 sta- 
tions and one near each end of the tuhe) and a self -balancing Indlcatlng- 
type potentiometer. Air total temperatures are measured hy thermocotq)les 
located In the central paissages of the IruLet and outlet mixing tanks 
of the test section as shown In figure 1. The Inlet mixing tank consists 
of two concentric passages so arranged that the air makes two passes 
through the t ank before entering the test section, while the outlet 
mixing tank Is comprised of four concentric; tanks. The passage of the 
air through the central chamber and concentric passages prevents radia- 
tion errors in the temperature measurements. Baffles are provided in 
the central passages to Insure thorough mixing of the air before Its 
temperature is measured. 
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Total-Pressure Measurements 

An opening for taUng total-pressure measurements across the tuhe 
is located 6 inches from the exit of the tuhe as shovn in figure 2. A 
hole having a 0.15-inch diameter throu^ "which a total-pressure prohe 
enters the tuhe "was drilled in the tuhe wall, at ri^t angles to "the 
static-presstire taps . A prohe actuator to move the prohe and measirre 
its distance in"to the tuhe "was fitted to a 6-inch length of tuhing at 
the opening. The location of the total-pressure prohe "With respect to 
the opening is sho"wn in figure 2(h) . The prohe tip has a flat opening 
"With a hei^t of 0.004 inch and a 0.002-inch vail and is made so that 
the tip Just clears the edge of the 0.15-inch hole in the test section. 
The total projected area of the prohe in the direction of flov vith the 

tip at the center of the tube is about l-g- percent of the area of the 

tuhe, hut the effective blocking area at the tip is considerably less 
because the maiTi pojrtlon of the prohe is do"wnstream of the tip. 


Total-Temperature Measurements 

The same opening in the tuhe used for taking total -pressure measure- 
ments is \ised for taking total -temperature measurements. Figure 2(c) 
shovs the location of the total-temperature prohe "with respect to the 
opening in the tuhe. One of the t"wo prongs of the prohe is made of 
chromel "while the other is made of alumel. The thermocouple of the prohe 
consists of 0.001-inch chromel and alumel vire hutt-velded between the 
prpngs "\diich were 0.07-inch apart. Temperatures as measured by the 
prohe were read on a self-balancing indicating-type potentiometer. 


METHODS 
Test Procedure 

Preliminary runs "with no air flow were first made to determine both 
heat loss through the test-section insulation and req.uired power settings 
for the end heaters. The tuhe was heated to various tenperatures by 
adjusting the electrical power input to the tuhe, and the power settings 
for the end heaters were adjusted to give a uniform wall temperature for 
each run. 

For each test run the flow control was adjusted to give "the desired 
flow rate and the electric-power control for the tuhe was set to give 
the desired tuhe-"wall temperature. The power settings to he "used for 
the end heaters at each tenpera"ture were determined from the preliminary 
runs . For the runs made the tube-wall temperature near the exit varied 
het"ween 300° and 1500° F", and the Reynolds ntunher "with fliiid properties 
based on the hulk temperature varied between 8000 and 500,000. The max- 
imum Mach number was about 0.6. 
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IKie foUoving quantities were measured for each run: air flow^ 

static pressvires at the -wall, outside tuhe-wall temperatures, inlet and 
outlet total hulk-air temperature, and electrical-power input to the 
tuhe. Total- pressure and temperature surveys at a section 6 inches 
from the exit of the tuhe "vdiere the flow was practically fully developed 
were taken during some of the tests . Measurements were taken at points 
between the wall in which the prohe opening was located and the center 
of the tuhe. Measurements made near the opposite wall were Inaccurate 
hecaAise of distTorhance due to the prohe. Total-pressure surveys were 
taken for tuhe-waH temperatures up to 1500° F, 'vdiereas total-temperature 
surveys were made only at low wall temperatures (300° F) because of the 
delicacy of the temperature prohe. 


Reduction of Experimental Data 

Local. heat-transfer rate . - The local heat -transfer rate per unit 
area from the tuhe to the air at a cross-section of the tuhe was calcu- 
lated from 



where the change of enthalpy AS was obtained from the hulk total tem- 
peratxires of the air at the inlet and outlet of the tuhe. The heat gen- 
eration per square foot of tuhe area G was obtained by measurement of 
electric current and tuhe resistance and the heat loss L was deter- 
mined by heating the tuhe to various temperatures with no air flowing. 

The quantities G, L, and S are all functions of tuhe temperature; the 
subscript loc means that the quantities were determined at the local 
tuhe temperatiore, whereas the subscript ave means that they were 
determined at the average tuhe temperatirre. The quantity in the first 
parenthesis is the average heat transfer per unit area; the quantity 
in the second parenthesis is a correction factor for local heat transfer. 

The local heat-transfer rate q^ was determined from equation (27) 
rather than from G^^qc “ ^loc tecaiise the measurements of inlet and out- 
let air temperature were, in general, more accurate than those of tuhe 
resistance and heat loss . The use of this equation also reduces errors 
due to possible noneq-ulllhrium conditions. The quantity in the second 
parenthesis in the equation did not vary greatly from 1 for most of the 
tests . 

Inside tuhe-waH and local hulk temperatures . - The inside tuhe-wall 
temperature tQ was obtained from the measured outside tuhe-wall tem- 
perature tg^ by the following equation derived in reference 4: 
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In equation (28) the assumptions are Eiade that heat is generated uniformly 
across the tuhe "wall thickness and that the heat flov at every point In 
the tuhe is radially Inward. For the present tests the radial ten^iera- 
ture drop through the tuhe wall was very small compared with the differ- 
ence between the inside wall and the air hulk temperature. 

The local total hulk-air temperature at a point was obta i ned by 
subtracting from the measured hulk temperature at the tuhe exit the tem- 
perature change produced by the heat added to the air in the portion of 
the tuhe between the exit and the point in question.. 

Physical properties of air . - In order to compare the experimental 
results with the analysis given in reference 2 , the viscosity and ther- 
mal conductivity were both assumed to he proportional to The 

specific heat and Prandtl number were assumed not to vary with tempera- 
ture. Substantially different results would be obtained if a different 
variation of thermal conductivity with temperatiare were used. A conduc- 
tivity proportional to t0»85 is sometimes giv^ in the literature, but 
the actual variation of conductivity i^ith temperature has not been exper- 
imentally determined at hi^ temperatures. 

Shear stress . - The shear stress at the wall for fully developed flow 
is related to the friction-pressure gradient by the equation 

"0 - - 1 (i)f, (25) 

The friction-pressure gradients were obtained by subtracting cal- 
culated momentimi-pressure gradients from the measured static-pressure 
gradients along the tube; the momentum-pressure gradients were calcu- 
lated from 


fe) ' (30) 

Equation (30) is exact only 'vdien the velocity profile is uniform; the 
effect of nonunlform velocity profile on calculated shear stress was 
checked for turbulent flow and was found to be negligible for the pres- 
ent tests. In equation (30 ) (>^ was found from the perfect gas law 


P = 


(31) 
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and dp^/dx was o^btalned "by differentiation of tlie perfect gas law to 
be 


^ 1 dp p 

dx ^t^ dx " gf^2 dx 


(32) 


In this eq.uatlon t^ was foimd from tbe following relation obtained from 
tbe equations of energy, continiilty, and state; 


% - 


- 1 +, 




~w^~R^^ 
gJ c^A p 


w^ 

X a2 2 

g J CpA p 


(33) 


This eq.uation can be more conveniently used for calculation purposes if 
tbe radical is expanded in a binominal series. Tbe q[uantity dt^/dx ^ras 
obtained by differentiating tbe expanded eqiiatlon for t^ id.tb respect 
to x; 


dx 


dx 



.2 r2 f 3pr/ ^ - 4T," g 


^ CpA 


where 


(34) 


dx wCp 


(35) 


Tbe static-pressure gradients in equations (32) and (34) were grapbically 
detennlned by plotting pressure against distance along tbe tube and 
drawing a tangent to tbe curve at tbe point in question. 


Velocity distributions . - For low air flow rates, incompressible 
flow theory was used; tbe velocities were calculated from tbe equation 
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P - p = -| pu^ (36) 

■where p \ras found from the perfect gas la'w p = pgRt and t 'VTas tahen 
equal to the total tenperatiire . In this and in all of the succeeding cal- 
culations, the static pressure "was assumed to he uniform across the tuhe. 

For Mach numbers greater than 0-2, velocities were calculated from 
the relation 


■^diere 


P 

P 


t 


r 



(37) 


For hoth lo'v/ and high Mach nunibers the total temperature in the 
eq_uations was obtained from the theoretical curves in figure '3. (Fig. 3 
is applied to compressible flow by replacing f*" by T^. This procedxire 
is legitimate for Prandtl numbers close to 1.) 


RESULTS AND DISCUSSION 

Predicted Temperatiire and Velocity Distributions 
(Prandtl Nuniber = 0.73) 

Temperature and velocity distributions for a Prandtl n\imber of 0.73 
as calculated from equations (lO), (n), (20), and (14) are plotted in 
figures 3 and 4. The values for the constants (n = 0.109 and x - 0.36), 
■which are fo\md in reference 3 from the experimental data for flow ■with- 
out heat transfer, are used for plotting the equations. In reference 3 
it is found that the equation derived for flow close to a wall ■^d.thout 
heat transfer agrees closely with the data for y+ less than 26 and that 
the equation for flow at a distance from a wall without heat transfer 
fits the data for y"*" greater than 26. For plotting the present equa- 
tions for flow ■\d.^t±L heat addition to the gas, the same limits of applic- 
ability for the equations for flow close to a wall and at a distance 
from a wall are used. It can be seen from the figures that the exact 
point of intersection of the curves representing the t^wo equations is 
not critical, especially for high values of 3, inasmuch as the slopes 
of the -two curves at their intersection do not differ greatly. 
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OOie exponent d in the equations vas found from Tiscosity data for 
air to have an average value of 0.68 for temperatiares between 0° and 
2000° F. Althou^ this value was obtained specifically for air, the 
values of d for most common gases do not vary greatly from this value^ 
so that the curves plotted shotiLd be applicable to most gases with 
Prandtl numbers close to 0.73. 

The plots of the equations in figures 3 and 4 Indicate that t'*' 
and u'*' are not equal as was the case for a Prandtl number of 1 (ref- 
erence 2) . Both profiles, however, show trends with increasing values 
of the heat-transfer parameter p •tJhich are similar to those found in 
reference 2 for a Prandtl number of 1. 


Nusselt Numbers 

In figures 5, 6, and 7, experimental and predicted Husselt numbers 
for air (Pr = 0.73) are plotted against Reynolds numbers for varlovis 
values of the heat-transfer parameter p. The experimental Nusselt 
numbers and Reynolds numbers were obtained at a point 6 inches from the 
exit of the tube idiere the flow was practically fully developed and 
were calculated from the definitions of Reynolds number, Nusselt number, 
and heat -transfer coefficient based on total-temperature difference, as 
given in the list of symbols. The predicted Nusselt and Reynolds num- 
bers were obtained from the equations given in the analysis for Prandtl 
numbers differing slightly from 1 and from the plots in figures 3 and 4. 

In figure 5 the physical properties, including density, in the 
Reynolds number and Nusselt number are evaluated at the static buht 
temperature. Both the experimental and predicted Nusselt numbers at 
a given Reynolds number show a decrease 'id.th increasing values of p 
or of ratio of wall to bulk temperature. In figure 6 the physical prop- 
erties, including density, in the Reynolds and Nusselt numbers are eval- 
uated at the wall temperature. For this case, both the experimental and 
predicted Nusselt nttmbers at a given Reynolds number increase with 
increasing p. In figure 7 the properties in the Reynolds and Nusselt 
numbers are evaluated at tQ which is sli^tly closer to the bulk 
temperature than the average of the wall and bulk temperatures . It is 
observed that the effect of p or of ratio of wall to bulk temperature 
on both the experimental and predicted Nusselt nimdbers is prsictically 
eliminated when the properties are evaluated at this temperature. The 
data follow the predicted line very closely for Reynolds numbers above 
15,000. For low Reynolds numbers the separation of the data from the 
predicted line is probably caused by a partial transition from turbulent 
to laminar heat transfer, idiich was not considered in the analysis. Ohe 
same trend at low Reynolds numbers was also observed in the data given 
in reference 1. It is seen that the transition region extends to higher 
Reynolds numbers for heat transfer than for friction. 
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It should he en^phasized that in oi^er for the foregoing correlation 
to holdj the same assumptions for the variation of physical properties 
id.th temperature must he made for calculating the ejcperimental results 
as were made in the analysis: constant specific heat and hoth thermal 

conductivity and viscosity proportional to If the conductivity 

had heen assumed proportional to t^*®^ as sometimes given in the lit- 
erature, the effects of ratio of wall to hilLk temperature on the Nusselt 
number correlation would have heen eliminated hy evaluating the fluid 
properties at a temperature close to the wall ten^perature . The actual 
variation of the thermal conductivity of air trith temperature has not 
heen experimentally determined at hi^ temperatures. In reference 1 
it is found that the hest correlation of average heat-transfer coeffi- 
cients for hoth heating and cooling was obtained hy assuming the thermal 
conductivity proportional to t*^*^ and evaluating the properties. 
Including density, in the Eeynolds and Kusselt mnnhers at the average of 
the wall and hulk teraperatvires . If a constant rather than a variable 
specific heat had heen used in reference 1, it would have heen necessary 
to use a conductivity proportional to tO-6 in order to eliminate the 
effect of ratio of wall to hulk temperature. The data for average heat- 
transfer coefficients given in reference 1 are therefore in substantial 
agreement with the present data and analysis for local heat-transfer 
coefficients. For the range of ratios of wall to hulk temperatures used 
in the present tests, as good a correlation could he obtained hy evalu- 
ating the properties at the average of the wall and hulk temperatures 
as hy evaluating them at to. 4* 

Although extensive data were not obtained in the low Reynolds num- 
ber range at high ratios of wall to hulk temperature, the data taken 
Indicate that in this range the effects of p are more nearly elimi- 
nated hy evaluating the properties at the hulk temperature than at 
tQ^^. This is in agreement with the analysis given in reference 5, 
where it is shown that for laminar flow the properties should he eval- 
uated at a temperature between the hulk temperatinre and the temperature 
at the center of the tube. It would he expected that in the transition 
region the properties should he evalviated at temperatures between this 
temperature and t^ 


Friction Factors 

In figure 8, experimental and predicted friction factors for air 
are plotted .against Reynolds h\imhers for various values of the heat- 
transfer parameter p. The physical properties, including density, 
in the friction factors and Reynolds numbers in the figure are evalu- 
ated at 4 ' The experimental friction factors and Reynolds numbers 
were obtained at a point 6 Inches from the exit of the tube and were 
calculated from the definitions of friction factor and Reynolds number 
given in the list of symbols . 
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The plot of predicted friction factors given in figure 8 indicates 
that the effects of p or of ratio of ■wall to hulk tenperature are 
nearly eliminated, hy evaluating the properties at tQ The data at 
high Reynolds numbers also Indicate that the ■ effects of ratio of wall 
to hulk temperature are eliminated hy evaluating the properties at that 
temperature. The data agree closely "with the predicted line in the 
high Reynolds number range. At lower Reynolds numbers the agreement is 
not as good, althou^ it may he ■within the experimental accuracy. Meas- 
urements of friction factors are, in general, much less accurate than 
those '’of heat-transfer coefficients because of difficulties in measuring 
static -pres sure gradients . In reference 1 it is found that some of the 
data for average friction factors at lo^w Reynolds numbers lie both belo^w 
and above the predicted line in figure 8, ■vdien the fluid properties are 
evaluated at the average of the wall and bulk temperatures . 


Velocity Distributions 

Experimental and predicted velocity distributions for various val-^ 
ues of p are sho^wn in figure 9. Rectangular rather than semllogarithmlc 
coordinates are used inasmuch as data are not shown for flow close to 
the wall; these data are not sho^wn because the distributions were meas- 
ured at hi^ Reynolds numbers where the severe velocity gradients and 
the presence of the hole in the tube wall make the accuracy of the meas- 
urements doubtful. 

It can be seen that the experimental and predicted values are in 
substantial agreement altho-ugh there is some scatter in the data. Inas- 
much as the experimental distributions were measiired in the high Reynolds 
nimiber range ■where the experimental friction factors fell on the pre- 
dicted line, the agreement between the experimental and predicted veloc- 
ity distributions might be expected. It appears that meas'urement of 
velocity distributions with heat transfer yields little Information in 
addition to that obtained by meeisurement of friction factors and heat- 
transfer coefficients other ■than possible Infoimatlon concerning changes 
in shape of the velocity profile ■with heat transfer. However, as indi- 
cated in figure 10, changes in profile shape should be very small over 
most of the tube cross section for the range of values of p obtained 
in the tests (O to 0.028) and probably could 3iot be measured 
experimentally . 


Temperature Distributions 

Experimental and predicted tenperature distributions are plotted 
in figure 11. Data for temperatxrre distributions close to the wall are 
sho^wn only at low flow rates because at high flow rates the severe 
velocity and temperature gradients and the presence of the hole in the 
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tube wall mate the accuracy of the measurements douhtful. The dlstri- 
hutlons shown for the region close to the wall may he subject to 
errors due to conduction along the thermocouple prohe prongs Inasmuch 
as the velocities In that region are very small . Temperature distri- 
butions were not measured for hl^ values of p or of ratio of wall 
to bulk temperature because of the delicacy of the temperature probe and 
because It appeared that such measurements would not yield appreciable 
Information beyond that obtained by measuring heat-transfer coefficients 
and friction factors at hl^ values of p. Changes In temperature pro- 
file shape with heat transfer should be small as In the case of velocity 
profile (fig. lO) . 

It Is seen In figure 11 that the experimental temperature distribu- 
tions agree fairly well with those predicted except In the low Reynolds 
number region where transition from turbulent to laminar heat transfer 
is taking place. This separation of data from the predicted line at 
low Reynolds numbers is not present in the velocity distribution data 
given in reference 3, so that it appears that the eddy dlffuslvities 
for heat arid momentum, transfer cannot be considered equal in the low 
Reynolds or Peclet niimber range. The eddy dlffusivlty for heat trans- 
fer in this range is probably a function of the thermal dlffusivlty 
V(pSCp) well as of the flow conditions. 


Predicted Effect of Variation of Shear Stress 
and Heat Transfer Across Tube 

A comparison between equations (7B) and (8B) (appendix B) ^ which 
take into aocount the variation in shear stress and heat transfer across 
the tube^ and the equations- from reference 2 ., iThlch assume uniform shear 
stress and heat transfer^ is shown in figure 12 for varioiis values of 
the heat-transfer parameter p. The intersection of the curves for flow 
close to and far from the wall is taken at y"^ = 30 rather than 26 for 
variable shear stress in order to make the mean deviation of equa- 
tions (7B) and (8B) from those given in reference 2 a mini-mum . It is 
seen that the effect of variable shear stress and heat transfer on the 
velocity distributions is sma.1 1 and is probably within the accuracy of 
flow and pressure measurements . The effect is greater for heat extrac- 
tion (negative p) than for heat addition to the gas . 

Althou^ the preceding calculations were carried out for a Prandtl 
number^ of 1, it is evident that the effects of- variable shear stress 
and heat transfer would also be small for a Prandtl number of 0.73. 
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Predicted Effect of Mclecular Shear Stress and Heat Transfer 
in Region at a Distance from Wall 

The results of the analysis of the effect of molecular shear stress 
and heat transfer in the region at a distance from the -wall are shown in 
figures 13 and 14. For the region close to the wall, the curves were 
calculated from the equation given in reference 2; for flow close to a 
wall for a Prandtl numher of 1 and for the region at a distance from the 
wall, ttie curves were calculated from equations (3C) and (4C) (appen- 
dix C) . The three experimental constants 'involved in the plots are n, 
X, and all of which are determined from experimental data for 

adiabatic turbulent flow from reference 3. The values for n and 
are taken directly from reference 3, and y^^^ is so determined that 
the curve for p = 0 agrees with the data for adiabatic flow. The 
slope v^^ in equation (3C) is set equal to the slope of the equation 
for flow close to the wall at (yq = 0.233 for p = 0.05). 

By including the moleciilar shear stress and heat transfer in the 
equation for flow at a distance from the wall, the region of applica- 
bility of the Karmm similarity expression is expended from y+>26 
to y"^>16. At y"^ = 16, the ciirve joins smoothly with the curve for 
flow close to the wall with no discontinuity in slope, that Is, the 
eddy dlffusivity as determined from the equations is continuoiis from 
the wall to the center of the tube. 

A conparison of the curves plotted with the molecular shear stress 
and heat transfer in the region at a distance from the wall neglected 
with those plotted with these factors considered in this region is pre- 
sented in figure 14. Only a small difference in results is obtained. 

As a result of this agreement, the heat-transfer coefficients and fric- 
tion factors will also be unaffected so that the molecular shear stress 
and heat transfer can be neglected with good approximation in the region 
at a distance from the wall. Although the preceding calculations were 
carried out for a Prandtl number of 1, it is evident that the same 
conclusions would hold for a Prandtl number of 0.73. 


SUMMARY OF ESBULTS 

The following results were obtained from the analytical and experi- 
mental investigation of fully developed turbulent flow and heat transfer 
in smooth tubes with variable fluid properties : 

li Substantial agreaiient was obtained between experimental and 
analytically predicted heat transfer and friction correlations eis well 
as velocity and tenperature distributions . • 
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2. 03ie test clieck of the analysis -was obtained hy measurement of 
local heat-transfer coefficients . Both the experimental and analytical 
results indicated that, for Reynolds numbers above 15,000, the effects 
of ratio of wall to bulk temperature on the Husselt number correlation 
can be eliminated by evaluating the fluid properties, including density, 
in the Reynolds and Russelt nimibers at a tenperature close to the aver- 
age of the wall and bulk tenperatures . In using the preceding result 
for calculations, it is important to use thermal conductivities and 
viscosities both proportional to the absolute temperature raised to the 
0.68 power and constant specific heat. 

3. The analysis Indicated that the effect of variation of shear 
stress and heat transfer across the tube on the velocity and tempera- 
ture distributions for turbulent flow with variable fluid properties 
is slight . 

4. The analysis indicated that the effect of molecular shear stress 
and heat transfer in the region at a distance from the wall is sli^t 
for tnrbulent flow with variable fluid properties . 


Lewis Fli^t Propulsion Laboratory 

National Advisoiy Committee for Aeronautics 
Cleveland, Ohio, October 18, 1951 
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APPENDIX A 
SYMBOLS 


A 


C, 


D 

d 

G 

g 


The folloviiig symbols are used in this report; 

cross-sectional area based on inside diameter of tube, sq. ft 
constants of integration 

specific heat of fluid at constant pressure, Btu/(lb)(*^) 
inside diameter of tube, ft 

exponent that describes variation of viscosity of fluid "with 
temperature 

heat generated in tube, Btu/(sec) (sq. ft) 
acceleration due to gravity, 32.2 ft/sec^ 


H 

h 


J 

k 

^b 

^0 

L 

n 

P 


enthalpy, Btu/lb 
heat-transfer coefficient. 


<30 


'^ 0 -% 


for compressible fluid or 


q_ u u 

r— for incompressible fluid, Btu/(sec) (sq. ft)(°F) 

■^0 “ 


mechanical equivalent *of heat, 778 ft-lb/Btu 

thermal conductivity of tube material, (Btu) (ft) /(sec) (sq ft) (°F) 

thermal conductivity of fluid, {Btu)(ft)/(pec)(sq ft)(°F) 

thermal conductivity of fluid evaluated at t^, 

(Btu)(ft)/(sec)(sq ft)(°F) 


thermal conductivity of flvild evaluated at irall, 
(Btu)(ft)/(sec)(sq ft)(°F) 

thermal conductivity of fluid evaluated at t^^ . , 
(Btu)(ft)/(sec)(sq ft)(^^) 

heat loss through insulation, Btu/(sec)(sq ft) 

constant 

total pressure, Ib/sq ft absolute 
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P 

q 

R 

^0 

3 

T 

Tb 

Tq or to 
t 

to. 4 
u 

% 

w 

X 

y 

r 

e 

c 

“h 

X 


static pressure, Ib/sq ft absolute 

rate of heat transfer toward tube center per unit eirea, 
Btu/(sec)(sq ft) 

rate of heat transfer at inside wall to-vrard tube center per 
unit area, Btu/(sec) (sq. ft) 

perfect gas constant, ft-lb/(lb ) (°R) 

inside tube radius, ft 

outside tube wall radius, ft ■ 

inside surface area, of tube , sq ft 

total temperature, °R 

buli or average total ten 5 >erature of fluid at cross section 
of tube, 

absolute wall temperatiire, 
absolute static temperature, °R 

bulk or average static temperatxire of fluid at cross section 
of tube, °R 

film temperature, 0.4(tQ -fb) 

V 

time average velocity parallel to axis of tube, ft/sec 

bulk or average velocity at cross section of tube , ft/sec 

fluid-flow rate, Ib/sec 

axi a1 distance from tube entrance, ft 

distance from tube wall, ft 

ratio of specific heats 

coefficient of eddy diffusivity for momentum, sq ft/sec 
coefficient of eddy diffusivity for heat, sq ft/sec 
Karman constant 
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1"0 

P 

PO 

PQ.4 

X 


absolute viscosity of fluids (it) (sec)/sq ft 
absolute viscosity of fluid evaluated, at t^^ (lb) (sec)/sq ft 
absolute viscosity of fluid at vail, (lb) (sec)/sq. ft 
absolute viscosity of fltiid evaluated at (lb)(sec)/sg ft 

mass density, Ib-sec^/ft^ 


2 / 4 

bulk or average density at cross section of tube, lb-sec /ft 
mass density of fluid at vail, Ib-sec^/ft^ \ 

density of fluid evaluated at t^ Ib-sec^/ft^ 
shear stress in fluid, Ib/sq ft 
shear stress in fluid at vail, Ib/sq ft 

Subscripts: 

fr on friction-pressure gradient' 

mom on momentiim-pressure gradient 

Dimensionless parameters: 

‘3o'\Ao/ PO 


^0 


P 


heat-transfer parameter. 


CpSToto 


0 


friction factor vith density evaluated at t^. 


D(dp/dx)p^ 2 t^ 


”0 


^0.4 


M 




2po^b P0% 

friction factor iri.th density evaluated at tg ^ 
Mach number. 


u 




Nusselt number vith thermal conductivity evaluated at t^^, 
hD 


WUq 


Nusselt number vith thermal conductivity evaluated at t^ 
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Nu, 


0.4 


Pr 


Nusselt nximber Tfith thermal conductivity evaluated at t 


0.4 


Prandtl numher. 


Cpl^g • 


Re^ 


Re 


0 


Re, 


'0 


0.4 

+ 




u+ 


Uv 


Un 


^1 


Reynolds number irlth density and viscosity evaluated at t^, 

i*b 

Reynolds number with density and viscosity evaluated at tQ 
Reynolds number with density and viscosity evaluated at tQ_^ 

aAo/^0 


tube -radius parameter. 


^^o/po 


static-temperatiore parameter. 


('faQ-'^)gpS^O 1- V~bQ 


‘3o'\Ao/ PO 
ll 

bulk static-temperature parameter, ■ 5 - 1 1 - x~~l 

p \ w 

value of t"*" at 


velocity parameter. 


u 




bulk; velocity parameter. 


value of u^ at yj_^ 


du 

wall-distance parameter. 


% 




^Ac/po 


value of y"^ at intersection of curves for flow close to 
wall and at a distance from wall 
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APPENDIX B 

ADTALYSIS INCLUDING EEFECT OF VARIATION OF SHEAR STRESS AND 
HEAT TRANSFER ACROSS TUBE (PRANDTL NUMBER = l) 

The following analysis wets male to determine the effect of varia- 
tion of shear stress and heat transfer across the tube on the velocity 
and temperature distributions with variable fluid properties . The var- 
iations of shear stress and heat transfer in the region close to the 
wall are negligible so that only the region at a distance from the wall 
is considered. 

The relation for the variation of shear stress with radius for fully 
developed flow is obtained by equating the shear forces to the pressure 
forces acting on a cylinder of fluid of arbitrary radius and differen- 
tial length. This relation gives 

The heat transfer per unit area varies in approximately the same 
way as the shear stress, elthou^ the variation is not linear. For the 
purpose of determining the general effect of the variation of heat 
transfer, the following relation is assumed: 

1 = qo (2B) 

Substituting equations (1B) and (2B) into equations (l) and (2), 
dividing equation (2) by equation (l), and integrating the result 
between the wall and a point in the fluid give, for a Prandtl number 
of 1, 


Cpg^O 



- t 


(3B) 


or 


= 1 - Pu+ (4B) 

With substitution of equations (4), (1B), and (4B), equation (l) 
can be written in dimensionless form with the molecular shear stress 
neglected as 
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( 1 y~^ \ = X (duVdy+) 

V -o7 (dV/ay+")2 


or 


d^u+ 


= - X 


dy-* 




/ du'* '^ 


(5B) 


tdiere tlie negative sign "was selected in taking tlie square root in order 
to make x positive. 

In order to solve equation (5B) let 


du"* 




<3y‘ 


Integration of equation (6B) gives 


= V 


(6B) 




y+ = y^+ + 


Un 


du~^ 

+ V 


(7B) 


With substitution of equations (6B) and (7B), eqization (5B) becomes 


dv 

V 


= - X 


du"" 


f(l- Pu+) 1 


1 I du+) 

V vJv V 


or 


-X 






du 


/(I - 1 - 


V = v^e 


yi _ 1 

^ ^ 


du~* 

IU]_+ V 


(8B) 


aZQ2. 
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Equation (8B) can "be solved iDy iteration by substituting aissumed values 
for and v into the ri^t side of the equation and calculating new 

values of y. After the relation between u'*' and v has been obtained 
for various values of P' and y* can-.be calculated from equa- 

tion (7B). 
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APPEKDIX C 


ANALYSIS HJCLUDING EFFECT CEF MOLECULAR SBEAR STRESS AND 
HEAT TRANSFER IN REGION AT A DISTANCE FROM WAIL 
(PRANDTL NUMBER = l) 

The effect of the molecular shear stress and heat-transfer terms 
in equations (l) and (2) on the velocity and ten 5 >erature distributions 
is investigated in this section. ' When these terms are retained^ the 
relation between temperatiore and velocity is given by equation (4B) as 
before. Equation (l) becomes, when written in dimensionless form for 
flow at a distance from the wall (Pr = l). 


1 = (1 


„ 4 _\d du'*' 
Pu+) ^ + 
dy^ 


x^(du+/dy+)^ 

2 

(l - Pu"*") (_d\‘^/dy''‘^) 


For solving this equation let 


du-*- 

dy-^ 


V 


(1C) 


(2C) 


Then 


d^u~*~ _ dv _ ^ dv 
dy+^ dy'*' du+ 

V 

As a result of substitution of these quantities in equation (ic) 


dv 

V 


V(l - (1 - pu+)'Vl(l - pu+) 

• i— mJ 


du' 


+ 


As in 
root. 


equation (l8) the negative sign ^Jas selected in taking the square 
Integration yields 



du+ 

/^[l - (1- pu+)'3i^(l- pu+) 


V 


(3C) 
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■where the suhscript 1 indicates that the values of and v are 

at the lo'west value of y"*" for Tdilch the equation applies. By assuming 
values of v coirespondlng to values of u+ and then calculating new 
values of v hy nvunerical ly integra’ting the right side of the equation, 
eq^uation (3C) can he solved hy Iteration. After the relation he-tween 
u"*" and V has heen obtained from eqtiation (3C), the relation he'tween 
u’*’ and y”^ can he obtained hy integrating equation (2C) as foUo-ws: 



■^•There y^"^ is the lowest value of y"*" for -vdiich equation (ic) applies. 
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Figure 8. - Variation of frlotlon factor with BoynoldB mniber for flow of air with heat addition an! properties evaluated 
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